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Abstract. We show that the Big Bang singularity of the Friedmann- 
Lemaitre-Robertson- Walker model does not raise major problems to 
General Relativity. We prove a theorem showing that the Einstein equa- 
tion can be written in a non-singular form, which allows the extension 
of the spacetime before the Big Bang. 

These results follow from our research on singular semi-Riemannian 
geometry and singular General Relativity [281 1291 131j (which we applied 
in previous articles to the black hole singularities [321 1331 l34l I30| ). 



0.1. The universe. According to the cosmological principle, our expand- 
ing universe, although it is so complex, can be considered at very large scale 
homogeneous and isotropic. This is why we can model the universe, at very 
large scale, by the solution proposed by A. Friedmann [HEIIH]. This exact so- 
lution to Einstein's equation, describing a homogeneous, isotropic universe, 
is in general called the Friedmann-Lemaitre- Robertson- Walker (FLRW) met- 
ric, due to the rediscovery and contributions made by Georges Lemaitre [19], 
H. P. Robertson [23 [Ml EZ] and A. G. Walker [SB]. 

The FLRW model shows that the universe should be, at a given moment 
of time, either in expansion, or in contraction. From Bubble's observations, 
we know that the universe is currently expanding. The FLRW model shows 
that, long time ago, there was a very high concentration of matter, which 
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exploded in what we call the Big Bang. Was the density of matter at the 
beginning of the universe so high that the Einstein's equation was singular 
at that moment? This question received an affirmative answer, under gen- 
eral hypotheses and considering General Relativity to be true, in Hawking's 
singularity theorem \1U\ [TT[ [T?j (which is an application of the reasoning of 
Penrose for the black hole singularities |22| . backwards in time to the past 
singularity of the Big Bang). 

Of course, given that the extreme conditions which were present at the 
Big Bang are very far from what our observations encountered so far, and 
our theories managed to extrapolate, we cannot know precisely what hap- 
pened then. If because some known or unknown quantum effect the energy 
condition from the hypothesis of the singularity theorem was not obeyed, 
the singularity might have been avoided, although it was a very high den- 
sity. One such possibility is explored in the loop quantum cosmology [Sj [3l [6] , 
which leads to a Big Bounce discrete model of the universe. 

We will not explore here the possibility that the Big Bang singularity is 
prevented to exist by quantum or other kind of effects, because we don't 
have the complete theory which is supposed to unify General Relativity and 
Quantum Theory. What we will do in the following is to push the limits of 
General Relativity to see what happens at the Big Bang singularity, in the 
context of the FLRW model. We will see that the singularities are not a 
problem, even if we don't modify General Relativity and we don't assume 
very repulsive forces which prevented the singularity. 

One tends in general to regard the singularities arising in General Rela- 
tivity as an irremediable problem which forces us to abandon this successful 
theory [13l [U [l6l [21 Hj. In fact, we will see that the singularities of the 
FLRW model are easy to understand and are not fatal to General Relativity. 
In [28] we presented an approach to extend the semi-Riemannian geometry 
to the case when the metric can become degenerate. In [29j we applied this 
theory to the warped products, and provided by this ways to construct ex- 
amples of singular semi-Riemannian manifolds of this type. We will develop 
here some ideas suggested in some of the examples presented there, and 
apply them to the singularities in the FLRW spacetime. We will see that 
the singularities of the FLRW metric are even simpler than the black hole 
singularities, which we discussed in [32 l [33l [3^ |3T]. 

0.2. The Priedmann-Lemaitre-Robertson- Walker model of the uni- 
verse. Let's consider the 3-space at any moment of time as being modeled, 
up to a scaling factor, by a three-dimensional Riemannian space (Sji^x;)- 
The time is represented as an interval I C M, with the natural metric — dt^. 
At each moment of time t £ I, the space T,t is obtaining by scaling (S, (7s) 
with a scaling factor a^(t). The scaling factor is therefore given by a func- 
tion a : / — )• M, named the warping function. The FLRW spacetime is the 
spacetime I x S endowed with the metric 

(1) ds2 = -dt2 + a2(t)dS2. 



BIG BANG SINGULARITY IN THE FLRW SPACETIME 



3 



It is the warped product between and (/, — dt^), with the warping 

function a : / — )• M. 

The typical space S can be any Riemannian manifold we may need for 
our cosmological model, but because of the homogeneity and isotropy con- 
ditions, it is in general taken to be, at least at large scale, one of the homo- 
geneous spaces S^, M'^, and H"^. In this case, the metric on S is, in spherical 
coordinates {r,9,(f)), 

(2) dT? = + r' {dO' + sin2 ed<i?) , 

where k = \ for the 3-sphere S^, A; = for the Euclidean space M'^, and 
k = —1 for the hyperbolic space 

0.3. The Friedman equations. Once we choose the 3-space S, the only 
unknown part of the FLRW metric is the function a(t). To determine it, 
we have to make some assumptions about the matter in the universe. In 
general it is assumed, for simplicity, that the universe is filled with a fluid 
with mass density p{t) and pressure density p{t). The density and the pres- 
sure are taken to depend on t only, because we assume the universe to be 
homogeneous and isotropic. The stress-energy tensor is 

(3) T-' = {p + p)u-u'+pg''\ 

where li" is the timelike vector field dt, normalized. 

From the energy density component of the Einstein equation, one can 
derive the Friedmann equation 

^a^ + k 

4 P= 

SttQ 

where k := — -r- (Q and c being the gravitational constant and the speed of 

light, which we will consider equal to 1 for now on, by an appropriate choice 
of measurement units). From the trace of the Einstein equation, we obtain 
the acceleration equation 

6 a 

5 p + 3p = . 

K a 

From these two equations we obtain the fluid equation, expressing the con- 
servation of mass-energy: 

(6) p = -3-{p + p). 

a 

Let's assume we know the function a. The Friedman equation Q shows 
that we can uniquely determine p from a. The acceleration equation deter- 
mines p from both a and p. Hence, the function a determines uniquely both 
p and p. 

From the recent observations on supernovae, we know that the expansion 
is accelerated, corresponding to the existence of a positive cosmological con- 
stant A |24^l23j. The Friedmann's equations were expressed here without A, 
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but this doesn't reduce the generahty, because the equations containing the 
cosmological constant are equivalent to those without it, by the substitution 



Therefore, for simphcity we wih continue to ignore A in the following, with- 
out any loss of generality. 




Figure 1 . The standard view is that the universe originated from 
a very dense state, probably a singularity, and expanded, with a 
short period of very high acceleration (the inflation). 

The standard view on cosmology today is that the universe started with 
the Big Bang, which is in general assumed to be singular, and then ex- 
panded, with a very short period of exponentially accelerated expansion, 
called inflation (Fig. [T]). 



The solution proposed here is simple: to show that the singularities of 
the FLRW model don't break the evolution equation, we show that the 
equations can be written in an equivalent form which avoids the infinities in 
a natural way. We consider useful to prepare the reader with some simple 
mathematical observations, which will clarify our proof. These observations 
can be easily understood, and combined they help us understanding the Big 
Bang singularity in the FLRW spacetime. 

1.1. Distance separation vs. topological separation. To understand 
the singularities in the FLRW model, it is useful to make a parallel with 
another type of singularities, and their standard resolution in mathematics. 
Let's consider a surface in M^. In general it can be defined locally as the 
image of a map / : [/ M^, where U £M? is an open subset of the plane. 
If the function / is not injective, the surface will have self-intersections. 
Another way to define the surface is implicitly, as the solution of an equation. 



(7) 



p —?■ p + K 
p — )• p — K 



lA 
lA 
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In this case it may happen again to have self-intersections. The typical 
example is the cone, defined as 

(8) -y^ -z^ = 0. 

We can desingularize it by making the transformation 

{X = u 
y =uv 
z = uw 

which maps the cylinder v'^ + w'^ = 1 from the space parametrized by 
(n, V, w), to the cone from equation Q, in the space parametrized by (x, y, z). 
This procedure is very used in mathematics, especially in algebraic geometry, 
and it is was studied starting with Isaac Newton [20]. 




Figure 2. The old method of resolution of singularities shows 
how we can "untie" the singularity of a cone and obtain a cylinder. 
This illustrates the idea that it is not necessary to assume that, 
at the Big Bang singularity, the entire space was a point, but only 
that the space metric was 0. 



The natural metric on the space (x, y, z) induces, by pull-back, a metric 
on the cylinder v"^ -\- w'^ = 1 from the space {u,v,w). The induced metric 
on the cylinder is singular: the distance between any pair of points of the 
circle determined by the equations u = and v'^ + w'^ = 1 is zero. But the 
points of that circle are distinct. 

From the viewpoint of the singularities in General Relativity, the main 
implication is that just because the distance between two points is 0, it 
doesn't mean that the two points coincide. We can see something similar 
even in Special Relativity: the 4-distance between two events separated by 
a lightlike interval is equal to 0, but those events may be distinct. 

1.2. Degenerate warped product and singularities. The mathemat- 
ics of General Relativity is a branch of differential geometry, called semi- 
Riemannian (or pseudo-Riemannian) geometry (see e.g. [21]). It is a gen- 
eralization of the Riemannian geometry, to the case when the metric tensor 
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is still non-degenerate, but its signature is not positive. In this geomet- 
ric framework are defined notions like contraction, Levi-Civita connection, 
covariant derivative, Riemann curvature, Ricci tensor, scalar curvature, Ein- 
stein tensor. These are the main ingredients of the theory of General Rela- 
tivity miEIlES]. 

The problem with the singularities is that there, these main ingredients 
can't be defined, or become infinite. The perfection of semi-Riemannian 
geometry is broken there, and by this, it is usually concluded that the same 
happens with General Relativity. 

In [28] we introduced a way to extend semi-Riemannian geometry to the 
degenerate case. There is a previous approach |171I18| . which works for met- 
ric of constant signature. Our need was to have a theory valid for variable 
signature (because the metric changes from being non-degenerate to being 
degenerate), and which in addition allows us to define the Riemann, Ricci 
and scalar curvatures in an invariant way, and something like the covari- 
ant derivative for the differential forms and tensor fields which are of use 
in General Relativity. After developing this theory, introduced in |28|, [3T] , 
we generalized the notion of warped product to the degenerate case, provid- 
ing by this a way to construct useful examples of singularities of this nice 
behaved kind [29j. 

From the mathematics of degenerate warped products it followed that 
a warped product like that involved in a FLRW metric (equation [T]) has 
only singularities which are well behaved, and which allow the extension 
of General Relativity to those points. At these singularities, the Riemann 
curvature tensor Rabcd is not degenerate, and it is smooth if a is smooth. 
The Einstein equation can be replaced by a densitized version, which allows 
the continuation to the singular points and avoids the infinities. 

1.3. What happens if the density becomes infinite? In the Friedmann 
equations Q, ([s]), and ([g]), the variables are a, the mass/energy density p and 
the pressure density p. When a — 0, p appears to tend to infinity, because a 
finite amount of matter occupies a volume equal to 0. Similarly, the pressure 
density p may become infinite. How can we rewrite the equations to avoid 
the infinities? As it will turn out, not only there is a solution to do this, 
but the quantities involved are actually the natural ones. As present in 
the equations, both p and p are scalar fields. But the adequate, invariant 
quantities actually involve the volume element, or the volume form 

(10) dvoi ■= ^/-gdt A dx A dy A dz, 

where by the factor ^/—g we mean ^/^^detg^. The densities are in fact not 
the scalars p and p, but the quantities pd^oi and pd^oi- They are differential 
4- forms on the spacetime, and the components of these forms in a coordinate 
system are p^/--g and respectively P\J—g- 

Another hint that the natural quantities are the densitized ones is given 
by the stress-energy tensor. When it is obtained from the Lagrangian, what 
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we actually get is the tensor density 

(11) T'^^^=_2^i%^ 

ogab 

The values p and p which appear in the Friedmann equations coincide 
with the components of the corresponding densities only in an orthonormal 
frame, where the determinant of the metric equals —1, and we can omit ^J—g. 
But when a — )• 0, an orthonormal frame would become singular, because 
det(7 — )• 0. A coordinate system in which the metric has the determinant 
— 1 will necessarily be singular when a(t) =0. In a non-singular coordinate 
system, det g has to be variable, as it is in the comoving coordinate system 
of the FLRW model. From ([T]) , the determinant of the metric in the FLRW 
coordinates is 

(12) det5 = -a^dets^s, 

where by dets gfs we denoted the determinant of the metric tensor of the 
3-dimensional typical space S. Since the typical space is the same for all 
moments of time dets (/s is constant. 

Given that the metric's determinant in the comoving coordinates is 

(13) V=5 = 

which tends to when a — )■ 0, we see that it might be possible for ^—g to 
cancel the singularity of p and p in p^^oi , respectively pd^oi ■ 

2. The Big Bang singularity resolution 



As explained in section ^1.3 we should account in the mass/energy density 



and the pressure density for the term \/—g. 

Consequently, we make the following substitution: 



P = PV-g = P'^^s/gr. 



(14) 

We have the following result: 

Theorem 2.1. If a is a smooth function, then the densities p, p, and the 
densitized stress-energy tensor Tab\/—g are smooth (and therefore nonsingu- 
lar), even at moments to when a(to) = 0. 

Proof. The Friedmann equation Q becomes 

(15) p= -a (d^ + k) ./g^, 

from which it follows that if a is a smooth function, p is smooth as well. 
The acceleration equation ([s]) becomes 

6 

(16) p + 3p = a^'dy/g^, 

K 

which shows that p is smooth too. Hence, for smooth o, both p and p are 
non-singular. 
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The four- velocity vector field is u = — , which is a smooth unit timelike 

Ot 

vector. The densitized stress-energy tensor becomes therefore 

(17) TabV^ = (p + P) UaUb + pgab, 

which is smooth, because p and p are smooth functions. □ 

Remark 2.2. We can write now a smooth densitized version of the Einstein 
Equation: 

(18) GabV^ + Agab\/^ = uTabV^. 

Remark 2.3. // a(0) = 0, the equation (15) tells us that p(0) = 0. From 
these and equation (16) we see thatp{0) = as well. Of course, this doesn't 
necessarily tell us that p or p are zero at t = 0, they may even be infinite. 
Figure\^ shows how the universe will look, in general. 



Figure 3. A schematic representation of a generic Big Bang sin- 
gularity, corresponding to a(0) = 0. The universe can be continued 
before the Big Bang without problems. 



Figure 4. A schematic representation of a Big Bang similar to 
an infinitesimal Big Bounce, corresponding to a(0) = 0, d(0) = 0, 
a(0) > 0. 
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Remark 2.4. One interesting possibility is that when a(0) = 0, also a(0) = 
0. In this case we may have a{t) > around t = 0, for example if d{0), and 
obtain a Big Bang represented schematically in Fig. This is very similar 
to a Big Bounce model, except that the singularity still appears. 

Remark 2.5. Note that we preferred the stress-energy tensor with lowered 
indices, Tab, io that with raised indices T"*, because the former involves the 
smooth metric gab, while the latter involves its inverse, g""^ , which is singular 
when gab becomes degenerate. Similarly for the Einstein Equation. The two 
versions are equivalent only when the metric is non- degenerate. 

3. The evolution of the universe 

This paper presented several scenarios concerning the Big Bang singular- 
ity, in the context of the Friedmann-Lemaitre-Robertson- Walker model. It 
was found that the singularity is of degenerate type, and the time evolution 
is not obstructed. The solutions are schematically represented in Figures [3] 
and|4j These models only tell what happens at the singularity. At a global 
scale, the universe may re-collapse in a similar singularity and then pass 
again beyond it, in a cyclic cosmological model, or may expand accelerating 
forever, as the present day observations seem to suggest [Ml [23] . Maybe the 
precedent universe, having t < 0, has no Big Bang at its origin, it just comes 
from the infinite past and collapses in a Big Crunch. Then, its Big Crunch 
becomes the Big Bang of our universe, and it starts its infinite expansion. 

The cosmological arrow of our universe points from the Big Bang toward 
the time direction where the universe expands, which is the direction in 
which t increases [t — t- -|-C!o). The cosmological arrow of the universe existing 
before the Big Bang points, of course, from the Big Bang, toward — oo (see 
Fig. §. 




Figure 5 . If the anterior universe has the cosmological arrow of 
time oriented toward our past, can we conclude that its entropic 
arrow of time also points toward our past? 

It is often assumed that the entropic arrow of time is determined by some 
special conditions existing at the Big Bang. Of course, the "entropic arrow" 
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may be undefined for a simple FLRW universe, bur it may be defined in 
universes which are at large scale approximated by FLRW models. If the 
cntropic arrow is determined by the cosmological arrow, then our model 
seems to suggest that the precedent universe has the entropic arrow of time 
oriented toward — oo, and its time flows from the Big Bang toward — oo, 
which is what the observers from the universe with t > would call past. 
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